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Abstract 

A universal exponential factor, 7 C = tt/2, is disclosed for the dimensional crossover of few-layer 
graphene (FLG) from two-dimensional graphene to three-dimensional graphite. "y c is found by 
analyzing available experimental data on different properties of FLG with varying thickness. A 
theoretical study on the phonon spectrum of the vertical acoustic mode in FLG is carried out to 
further check this exponential factor ~f c . Interestingly, the same exponential factor appears in the 
dimensional crossover of the phonon mode. It turns out that the exponential factor 7 C is related 
to the homogeneous Helmholtz-like molal equation in the mass transfer with a first order chemical 
reaction. The finding should provide valuable information for experimentalists and theorists in the 
future investigation on thickness dependent properties of FLG. 
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I. INTRODUCTION 



The advancement of producing few-layer graphene (FLG) of high quality enables FLG 
to be a very promising candidate for future nano-devices. The single layer graphene sheet 
(N = 1) is an one-atom-thick two-dimensional (2D) system; while FLG with N — > oo is a 
three-dimensional (3D) graphite system. Hence, the FLG provides a ready-made platform 
for the investigation of dimensional crossover from 2D to 3D by increasing the layer number 
N. It is quite interesting to investigate the evolution for different physical properties of 
FLG during the dimensional crossover from 2D graphene to 3D graphite. In a very recent 
experiment, [lj the in-plane thermal conductivity was measured for FLG with different N, 
and it was found that the dimensional crossover from 2D to 3D is very fast. For FLG with 
N > 4, the value of the thermal conductivity is already very close to that of the 3D graphite. 
There are some other available experiments on the dimensional crossover of different aspects 

n 

for FLG. Those experiments include the interface thermal conductance, |2| the Raman G 
mode,^ the second order D peak, [4] the surface roughness, Q the electronic band gap, j^l 
the surface potential, [3] etc. The main purpose of this work is to extract the common rules 
in the dimensional crossover of FLG from existing experiments, check them theoretically, 
and understand the underlying mechanism. 

In this work, by analyzing various experimental data, we find a universal exponential 
factor 7 C = 7r/2 for the dimensional crossover of FLG from 2D graphene to 3D graphite. We 
then perform theoretical check on 7 C by studying the phonon spectrum of the vertical acoustic 
(ZA) mode in FLG. We find that the phonon spectrum of the ZA mode varies exponentially 
during the dimensional crossover from 2D to 3D with exactly the same exponential factor 
7 C . Finally, we give a physical interpretation for the value of 7 C by relating it to the mass 
transfer phenomenon with a first order chemical reaction. 



II. ANALYSIS FOR EXPERIMENTAL RESULTS 



In this section, we analyze several experiments on the dimensional crossover from 2D 
graphene to 3D graphite for different physical properties of FLG. These experiments cover 
thermal properties, electronic properties, Raman studies, and mechanical property of the 
FLG. We will demonstrate that all experimental data can be well described by exponential 
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functions with a common exponential factor 7 C = ir/2. 

Fig. [I] (a) shows the in-plane thermal conductivity of the FLG with different layer number 
N. Both the value and error bar in the figure are from Fig. 2 (b) in Ref. Q. The temperature 
dependence of the G mode was used to detect the temperature gradient in the FLG. The 
obtained temperature gradient together with the thermal power of the laser can be used 
to derive the value of the thermal conductivity. This method has been applied in single 
layer graphene sheet by the same group and a thermal conductivity of 4000 W/m/K was 
discovered in this 2D system. [8] The thermal conductivity of the regular graphite is about 950 
W/m/K. Fig. [T] (a) exhibits the dimensional crossover of the thermal conductivity from the 
2D graphene to the 3D graphite systems. The thick FLG sample with N = 8 has a larger 
uncertainty in its layer number with 6 < N < 10. The measured thermal conductivity 
is lower in thicker FLG samples. This result was attributed to the enhancement of the 
three phonon scattering channels by the inter-layer van de Waals (VDW) interaction. The 
experimental data can be fitted to an exponential function k = 992.2 + 5916.9e~ Ar / 7c with 
the exponential factor 7 C = tt/2. In the limit of large N, k — 992.2 W/m/K yields the 
thermal conductivity of the 3D graphite within five percentage. The small exponential factor 
7 C = tt/2 indicates a fast crossover from 2D graphene to 3D graphite with the increase of 
N. It means that the thermal transport capability of FLG with iV = 2 is already very close 
to 3D graphite. In following, we will show that the same exponential factor 7 C = tt/2 exists 
in other experiments. 

Fig. [1] (b) is the vertical thermal conductance per unit area (G) v.s layer number N. 
These data are for 3.8 fim FLG samples from Fig.3 (b) in Ref. 2l The thermal conductance 
was measured by time-domain thermoreflectance method, which may give a 20% uncertainty 
in measured value. For each N, there are two values of G corresponding to different exper- 
imental conditions. The low thermal conductance in the experiment was attributed to the 
Kapitza thermal resistance of various interfaces in this sandwich structure. Fig. [T] (b) shows 
that the thermal conductance decreases exponentially and reaches a saturate value quickly 
with the increase of layer number N. The data can be fitted to an exponential function, 
yielding the common exponential factor 7 C . This small exponential factor indicates that a 
FLG sample with N = 2 can be used to mimic the 3D graphite quite reasonably. 

Fig. [1] (c) shows the Raman study of the G mode in FLG with different N, which is from 
Fig.3 and Fig. 4 of Ref. [j. The error bar (about 0.01 cm -1 ) is an estimation of possible errors 
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introduced when data are read out from the experimental work. For each N, the variation 
in experimental frequencies is about 0.5 cm -1 , which was attributed to the spot-to-spot 
variation within the FLG films. [3 1 ] Fig. [T] (c) shows a red-shift of the G mode with increasing 
N. This red-shift is originating from the polarization of the bond charge due to the inter- 
layer VDW interaction of the FLG. [9] The experimental data in the figure are fitted to an 
exponential function. We find that the common exponential factor 7 C also presents itself in 
this experiment. In large N limit, the exponential function gives 1582.04 cm -1 for G mode. 
This value meets the frequency of G mode in 3D graphite (1581.91 cm -1 ). This agreement 
confirms that the exponential function with common factor 7 C is a good description of the 
experimental frequency of G mode. 

The position of the second-order Raman D peak in FLG is shown in Fig. [T] (d). The data 
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are read out from Fig. 2 of Ref. U. Similar experimental results can also be found in Ref. 
This frequency corresponds to the second-order Raman scattering of zone-boundary phonon 
around 2700 cm -1 . For the second-order Raman D peak, there is only a single component 
in graphene, while two components are found in graphite. J4] We have extracted the position 
of this Raman peak by reading the value of the frequency at the center of the Raman peak. 
In the experiment, two Renishaw spectrometer were applied: 514 and 633 nm. We focus on 
the experimental data under 514 nm. Compared with results at 514 nm, the shape and the 
position of the Raman spectra at 633 nm show an overall shift, but the iV-dependence is 
similar. As the Raman spectrum experiment is of high accuracy, the error bar (about 0.7 
cm -1 ) for each data in Fig. [1] (d) is estimated from possible errors during we reading out 



data from Ref. 
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This figure shows that the frequency of the 2D mode is well described 
by to = 2726.4 - 79.6e _iV/7c , where the constant 2726.4 cm 1 is exactly the value in 3D 
graphite. Interestingly enough, the dimensional crossover of the frequency is determined by 
the common exponential factor 7 C = n/2, the same as previous experiments. 

As shown in Fig. [2] (a) is the roughness of FLG with different thickness. This experiment 
was performed by Lauffer et al. in Ref. 5). The interface-induced roughness R reflects the 
height variations of the FLG. The error bars in the figure are mainly caused by the variations 
in the tip geometry. The single layer graphene sheet has a rather large roughness _R=30pm, 
and R decreases quickly with the increase of N. Similar roughness of single layer graphene 



sheet was also observed by high-resolution atomic force microscopy. 
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The decrease of the 



roughness with increasing thickness was also found by Brar et al. [12], where the Bernal 
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stacking structure of FLG was observed. This experiment demonstrates that thicker FLG 
samples are more smooth due to the inter-layer VDW interactions. The roughness R can be 
fitted to exponential function R = 55.le~ N ^ c . This function yields a zero roughness for 3D 
graphite. Again, the common exponential factor 7 C exists in this experiment. 

Fig.[2](b) shows the electronic band gap v.s N in FLG. The angle-resolved photo emission 
spectroscopy was taken for FLG to deduce the electronic band gap in Ref. |6j. A band gap of 
0.26 eV is obtained for single layer graphene sheet, which benefits its application in electronic 
devices field. The value of E g = 0.025eV for N = 10 is actually the band gap in 3D graphite 
which was estimated from a band-structure calculation. 16|, ll7| This value should be set as 
the band gap for infinite large N. However, we find that N = 10 is already large enough. 
There is no difference in the fitting function if we set this value as the band gap for FLG with 
N = 100. As shown in the figure, the band gap decreases exponentially with the increase 
of layer number N. The common factor 7 C appears in this experiment as the exponential 
factor. Although the existence of this band gap is widely accepted, the origin of the gap is 
a current object of an ongoing debate. Zhou et al. interpreted this band gap as resulting 
from the breaking of sublattice symmetry due to the graphene substrate interaction in the 
epitaxially grown FLG samples on SiC substrate; [13j while Rotenberg et al. believed that 
the origin of the band gap is due to the modulation of the lateral structure of the FLG 
samples, [l^] A more definitive interpretation of the opening has been presented by Bostwick 
et al., where a renormalization of the bands around the Dirac crossing is caused by the 
electron-plasmon coupling, [jj]] 

The surface potential of FLG with different N is displayed in Fig. [2] (c). These experi- 
mental data are reported in Fig. 3 of Ref. ?j- Similar results can also be found from Ref. 
The surface potential is the electrostatic potential measured at the exterior layer of the FLG 
sample. The value was probed by the electrostatic force microscopy. It was proposed that 
the surface potential is affected by the charge exchange between the FLG and the substrate. 
We can see that the evolution of the surface potential with N can be well described by the 
exponential function, and we obtain the common exponential factor 7 C . 
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III. DIMENSIONAL CROSSOVER FOR ZA PHONON MODE 



In the above, a universal exponential factor 7 C = 7r/2 is revealed by analyzing various 
experimental results on the dimensional crossover from 2D graphene to 3D graphite. The 
rest of this paper is devoted to the discussion of the dimensional crossover for the ZA mode. 
ZA mode is the acoustic phonon in the vertical direction of the graphene. This is the so- 
called flexure mode in graphene which has a quadratic phonon spectrum due to the 2D 



plane sheet structure of graphene. 
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20| We study the evolution of the ZA mode during 



the dimensional crossover from 2D graphene to 3D graphite. The phonon dispersion of ZA 
mode can be described by a formula u = ak 13 , where a and (3 are two parameters. The 
phonon dispersion is obtained by solving the lattice dynamics of FLG at the equilibrium 
structure from the "General Utility Lattice Program" (GULP) 21]. The in-plane carbon- 
carbon interaction is described by the Brenner potential. 22] For the inter-layer interaction, 
since the distance between two layers is out of the interaction range of Brenner potential, 
we introduce Lennard- Jones potential, V(r) = 4e((a/r) 12 — (cr/r) 6 ), with e = 2.5 meV and 
a = 3.37 A. The value of the length parameter a is fitted to the space between adjacent 



layers in three-dimensional graphite 



23j as 3.35 A. The energy parameter e is fitted to the 



phonon dispersion along TA direction in graphite. 24J The cutoff for Lennard- Jones potential 
is chosen as 20 A, which is large enough; i.e., all results do not vary with further increasing 
this cutoff value. Fig. [3] (a) is the ZA mode for N = 1, i.e single layer graphene sheet. 
The ZA mode in graphene is a perfect flexure mode with (3 = 2. We have considered the 
long wave phonon with the wave vector from V to 0.03FA" in the Broullin zone, where 
TK = 47r/(3y / 3&) and the bond length b = 1.42 A. Panel (b) is the ZA mode in 3D graphite, 
where a = 1.8 and (3 = 1.0. Panels (c) and (d) display the two parameters a and (3 as 
a function of layer number 1 < N < 10. These two curves can be fitted to exponential 
functions, and the exponential factor 7 C = tt/2 is the same as what we have discovered from 
experiment analysis. In the limit of large N, a = 1.72 and (3 = 1.07 give the results for 3D 
graphite within some small errors. 

We point out that the difference in ZA mode for different N originates from the two 
nonequivalent sublattices in the honeycomb structure of graphene. In the elastic model, the 



FLG are constructed by many uniform elastic plates coupled through VDW interactions. 25] 
The solution for the ZA mode is u n (x,y) = Ae^ kxX+kyy ^ e wt ; which is independent of the 
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layer index n. This solution simply indicates that all graphene layers vibrate with the same 
amplitude and phase. As a result, the VDW has no contribution to the ZA mode, and the 
phonon dispersion of the ZA mode in FLG should be the same as single layer graphene, 
i.e to = 0.7 * k 2 cm -1 . However, the elastic model loses all lattice structure information, 
especially the nonequivalent property of the two lattice sites A and B. In the Bernal- 
stacked FLG, two adjacent layers have a relative shift by a carbon-carbon bond length. This 
relative shift can be considered by modifying the solution of the ZA mode to be n-dependent: 
U2n(x,y) = Ae l( - kxX+ky ^ ] e luJt for even layer and u 2n +i{x, y) = Ae % ^ kxX+kyV) e wi e %k ' h = u 2n e %k ' h for 
odd layer. The vector b is the relative shift between two adjacent layers. This solution 
exhibits that the two adjacent layers have a phase difference during vibration, which has 
effect on the long wave length phonon. As a result, the VDW interaction can affect the 
phonon dispersion of the ZA mode; thus the ZA mode is different in FLG with different N. 



IV. UNDERLYING MECHANISM 



The above universal exponential factor 7 C = ir/2 is probably rooting in the stacking con- 
figuration structure of the FLG, as this constant exists in quite different physical quantities. 
Let us consider a general quantity Q siting on one carbon atom in the first layer of the FLG. 
There is a possibility for Q to hop to other sites. We focus on the hoping phenomenon in 
the vertical (z) direction of FLG. Under the first-nearest-neighboring tight-binding approx- 
imation, the hopping parameter K in z direction is inverse proportional to the inter-layer 



space c. 



261 ] In the meantime, larger in-plane bond length b will result in larger K , because 



in-plane hoping is suppressed in case of larger b which enhances hoping in z direction. These 
two effects determine the hoping parameter as K = b/c. As a result, the hopping rate in 
z direction is g'q = —KqQ under the lowest-order approximation. Before writing down the 
equation for Q, it is important to point out that the hopping process discussed here is actu- 



ally analogous to the mass transfer phenomenon with a first order chemical reaction. |27| The 
quantity Q corresponds to the distribution of molecule concentration. Hoping in z direction 
is equivalent to the chemical reaction to absorb the molecule. The hoping rate g'q = —K Q Q 
is corresponding to the absorption rate of molecule under the first-order chemical reaction 
with Kq as the reaction constant. By analogy with mass transfer, the distribution of quantity 
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Q along z axis in steady state is governed by the molal equation: 27] 

d 2 Q 



dN 2 



- K Q = 0, (1) 



where the coordinate z is discretized in FLG system. The solution for this homoge 



neous Helmholtz-like equation under vanishing boundary condition is:|28| Q oc e^ N ^ c with 
7 C = yJ\/K$. This exponential solution is the origin of previous exponential functions 
in the dimensional crossover of different physical quantities; and the exponential factor 
7 C = y/l/Ko = a/c/6 = 1.54 agrees quite well with above 7 C = ir/2 = 1.57 from exper- 
imental and theoretical analysis. It should be noted that the we are studying the Bernal 
stacking FLG sample in our work. For turbostratic graphite or epitaxial graphene grown 
on SiC(OOO-l) with a peculiar turbostratic, non Bernal stacking, many of its properties are 
strongly affected. For example, the Raman G modes undergo stronger red-shift in AAAA 
stacking FLG than that of the Bernal stacking FLG;|9( thus the exponential factor is slightly 
smaller than tt/2. In our model, we have obtained the hoping parameter within the first 
order approximation. This approximation is more reasonable for Bernal stacking structure, 
where only one of the two sublattice has inter-layer neighbors. 

V. CONCLUSIONS 

To summarize, we find a universal exponential factor 7 C = 7r/2 during the dimensional 
crossover of FLG from 2D graphene to 3D graphite. We obtain this exponential factor 
through the analysis of several experiments on the transition of different physical properties 
of the FLG with increasing layer number N. We then study theoretically the ZA mode in 
FLG with phonon spectrum u = ak^ . We find that the two coefficients, a and /3, vary 
exponentially with the increase of N, and the exponential factor is exactly the same as what 
we have obtained from the analysis of experiments. We also demonstrate that 7 C is related 
to the mass transfer phenomenon with a first order chemical reaction. 
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FIG. 1: Experimental results on various physical quantities in FLG with different N. The constant 
7 C = 7r/2 is the common exponential factor, (a) In-plane thermal conductivity k measured in Ref.Q. 
(b) Interface thermal conductance G from Ref. 0. (c) The Raman G mode from Ref. 0. (d) The 
second order D peak in Ref. [4. 
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FIG. 2: Experimental results for the dimensional crossover from 2D to 3D of three physical quanti- 
ties in FLG. The constant j c = tt/2 is the common exponential factor, (a). The surface roughness 

). (| (c)^ The surface potential of FLG. [3] 
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FIG. 3: Phonon dispersion of ZA mode in FLG is ui = ak 13 . (a). ZA mode in graphene sheet is 
flexure mode with parabolic spectrum, j3 = 2.0. (b). ZA mode in 3D graphite is linear, /3 = 1. 
(c). The coefficient a v.s N. (d). The power factor j3 v.s iV. 
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